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Abstract: We study general two-body decays of primitive and non-primitive |-BPS 
dyons in four- dimensional type IIB string compactifications. We find a "master equa- 
\ tion" for marginal stability that generalises the curve found by Sen for |-BPS decay, 

^ and analyse this equation in a variety of cases including decays to |-BPS products. For 

i-BPS decays, an interesting and useful relation is exhibited between walls of marginal 
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stability and the mathematics of Farey sequences and Ford circles. We exhibit an 



example in which two curves of marginal stability intersect in the interior of moduli 
space. 
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1. Introduction 

In the last couple of years there has been renewed interest in the properties of dyonic 
black holes in four dimensions, particularly those associated to M = 4 compactifications 
(type II strings on K3 x T 2 or heterotic/type I strings on T 6 , as well as supersymmetry- 
preserving orbifolds of these systems) 0, g, |, |, g, g, Q, §, |0|, 0, [14]. 



A key advance has been a better understanding of a classic degeneracy formula due 
to Dijkgraaf, Verlinde and Verlinde[]T3[. Among other things, the generalisation of 
this formula to CHL orbifolds and the origin of a genus-2 modular form have been 
illuminated in many of these works. 
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Recent work has focused on the issue of marginal stability of these dyons. Curves of 
marginal stability for specific decays have been obtained!!, the impact of such decays 



on the degeneracy formula has been studied || [10], [13] and the decays across such 
walls have been identified with the disappearance of two-centred black holes from the 
spectrum [|12|, jl3|, following previous ideas in the M = 2 context [IB] . A formula has 



been proposed in [|13j to count the "immortal" dyons which exist everywhere. And very 
recently, Sen has considered the case of primitive dyons decaying into |-BPS states |H 



and demonstrated that this takes place only on surfaces of codimension 2 in moduli 
space. 

Clearly there is much more to be learned about this system. Among various inter- 
esting questions is a complete understanding of the possible marginal decays of |-BPS 
dyons, the impact of such decays on the degeneracy counting function, the role of multi- 
centred black holes in the decays, and the relevance of "non-primitive" dyons (which 
are related to Riemann surfaces of genus g > 2) to the counting problem. 

In the present work we take a step towards resolving the first problem. We consider 
the most general decay of a ^-BPS dyon into two decay products, each one of which can 
be either |- or ^-BPS. We find a necessary condition for marginal two-body decays and 
study the resulting equation in a variety of cases. It turns out that some solutions of our 
equation are "spurious" in the sense that they describe an inverse decay process rather 
than the forward decay. This puts constraints on the possible decay products which 



are identical to those found in |TJ] . We are also able to reproduce some of the results 
in Refs. IBIl ClS db special well as generalise them to the case of "non-primitive" 

dyons. On the way we will see that a known mathematical construction, that of Farey 
sequences and Ford circles, bears a remarkably close relation to the circles of marginal 
stability in Ref. and helps us understand the properties of these circles. 



2. The system 

We consider type IIB string theory compactified on K3 x T 2 . The resulting four- 
dimensional system has 28 U(l) gauge fields and their electric-magnetic duals. There- 
fore we can have dyons of charge (Q, P) where the first entry is a 28-component vector 
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denoting electric charge under these gauge fields and the second denotes the magnetic 
charge. The dyons will be |-BPS if the vectors Q, P are parallel, and ^-BPS otherwise. 

Note that a modular transformation of the 2-torus T 2 that changes its modular 
parameter as: 

r ?Z±* (2.1) 
cr + d 

with 

''' ^ e PSL(2, Z) (2.2) 



(2.3) 




sends the dyon charges to: 

'aQ + bP s 
K cQ + dP; 

We are interested in the marginal decays of these |-BPS dyons. The stability or 
otherwise is dictated by the charges carried by the dyons as well as the values of the 
moduli of K3 x T 2 . These are encoded as follows. Define the matrix: 

L = diag(l 6 ; (-1) 22 ) (2.4) 

In 4 dimensions there are, first of all, 132 moduli that can be assembled into a matrix 
M that is symmetric and orthogonal with respect to the L metric: 

M T = M, M T LM = L (2.5) 

The relevant inner product for an electric charge vector, which we will call Q 2 or Q ■ Q, 
is 1 : 

Q 2 = Q T (M + L)Q (2.6) 

Correspondingly we have: 

P 2 = P T (M + L)P 
P-Q = P T (M + L)Q (2.7) 

We will also make use of the quantities Qr, Pr defined such that 

Q R = Q T r Qr = Q T (M + L)Q (2.8) 



1 Because our focus is on microstates, our inner products are always dcfinccl with respect to the 
moduli at infinity, so this notation should not cause confusion. 
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and similarly for the other inner products (for details see for example [fL3|, [14]). 

In addition to the moduli appearing in M, there is the modular parameter of the 
4-5 torus: 

r = n+ ir 2 (2.9) 
The BPS mass formula for general |-BPS dyons is [|I7], 



Mbps(Q,P) 2 = -^(Q -fP) ■[()- rP) - 2 s /7-> X /MQ. P) (2.10) 



2 



where 

A(Q,P) = Q 2 P 2 -(P-Q) 2 (2.11) 

Before going on, it is useful to transform the dyon charges to bring them into a 
standard form. Consider the electric and magnetic charge vectors Q, P of the dyon and 
define [0]: 

/(<§, P) = gcd(Q A P) = gcd(Q i P j - Q j P l ), all i, j (2.12) 

If for a given dyon we find that I(Q,P) > 1, we first perform an SL(2,Z) transfor- 
mation as in Eq. fl2.3| ). Using some properties of finitely generated algebras (see for 
example Ref.[|19|, Chapter I, Section 8), we can always find such a transformation 2 that 
yields new dyon charges of the form (mQ' , nP 1 ) for some positive integers m, n and some 
new vectors Q' , P' such that I(Q', P') = 1. Under this transformation I(Q, P) remains 
invariant, so m, n must be such that I(Q, P) = mn. If the m, n so obtained are not 
co-prime then the dyon with those m, n will be marginally unstable at all points of 
moduli space. This does not mean a bound state does not exist, but that determining 
its existence is more subtle and requires actually quantising the system. Therefore we 
will restrict ourselves to the case where m, n are co-prime. 

To summarise, in what follows we assume that our dyons have charge vectors 
(mQ,nP) with co-prime m,n and with I(Q,P) = 1. The special case (m,n) = (1,1) 
will be called a primitive dyon. 



2 We are grateful to Nitin Nitsure for helpful discussions on this point. 
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3. Decays into a pair of dyons 

We can now examine the decay of a ^-BPS dyons into two other dyons. From charge 
conservation, the most general decay is of the form: 





(3.1) 

where Qi,Pi are arbitrary vectors in the (6, 22)-dimensional integral charge lattice. 



From the study of BPS string junctions and networks p0|, p2f , we know that the 
decay products can be mutually BPS with each other and with the initial state only if 
the corresponding charges all lie in a plane rather than being generic 28-dimensional 
vectors as above. However, the properties of the networks are determined in the present 
context not by the charge vectors Q, P but by their projections Q R , P R . Indeed it is 



only the latter which appear in the BPS mass formula Eq. ( [2.10| ) that we will be using. 



Therefore the BPS condition requires that the R projections of the final-state charges 
are in the same plane as those of the initial-state charges. This happens automatically 
in some cases, while in others it requires adjusting the moduli in M to make this 
happen. 

It follows that we must have the relation: 




miQ R + riP R \ + / m 2 Q R + r 2 P R \ ^ 
siQ R + jt-iPr J \s 2 Q R + n 2 P R J 



where the coefficients m^, n^, r^, Si satisfy: 

m 1 + m 2 = m, ni + n 2 = n, r% + r 2 = s% + s 2 = (3.3) 

We cannot, however, assume that these coefficients are integers since the above equation 
refers not to the original vectors in the integral lattice but to their projections to the 
Q R , P R plane. 

Without any additional conditions on these coefficients the decay products will 
both be |-BPS dyons. It is possible to have one or both of them be |-BPS by suitably 
constraining the integers, as we will see shortly. 

If M, Mi, M 2 denote the BPS masses of the initial state and the two decay products 
(for simplicity we henceforth drop the subscript BPS), we can use Eqs.( |2.10| ) and 
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(|3.2|) to evaluate the condition on the moduli imposed by the marginality condition 
M = Mi + M 2 . Because of the square root in Eq. ( |2.10| ), this is most easily done 
by computing a combination of squared masses that vanishes when the marginality 
condition is satisfied. 

First, define the angles 9 and 9i 2 by: 

9 = tan" 1 - 

Q -1 Qr ■ Pr fn a\ 

9 QP = cos (3.4) 

where Qr = \Qr\, Pr = \Pr\- Geometrically, 9 is the opening angle of the torus while 
9qp is the angle between the projected electric and magnetic charge vectors (which 
coincides with the angle appearing in the string junction description of the dyon). We 
also define a "cross-product" between the integers mi, n-i, m 2 , n 2 as: 

m An = min 2 — m 2 ni (3-5) 

Let us now find the condition that, at some point (s) in moduli space, the decay 
Eq. ( |3.2|) becomes marginal: M = Mi + M 2 . The BPS formula Eq. Q2.10 ) involves a 



square root on the RHS and another square root to extract M from M 2 . The simplest 
square-root-free expression that vanishes when the marginality condition is satisfied is 
the combination: 



M 4 + M 4 + M:; - 2{M 2 M\ + M 2 M 2 + M 2 M. 



) 



= (M- Mi- M 2 )(M + Mi + M 2 )(M - Mi + M 2 )(M + Mi - M 2 ) (3.6) 

Now we require this expression to vanish. However, subsequently we must check that 
on the vanishing curve, it is really the first factor of the RHS of Eq. ( |3.6j ) that vanishes 
rather than any of the other factors. Notice that the second factor never vanishes (since 
all the M's are positive), while vanishing of the third or fourth factor corresponds to the 
inverse decays Mi = M + M 2 and M 2 = M + M x . When we turn to a detailed analysis 
of marginal decay processes, it will be necessary to rule out these inverse decays before 
concluding that we are dealing with the correct decay mode. Only in the case where 
both the final products are |-BPS, this check becomes unnecessary because the reverse 
process is forbidden: a |-BPS state cannot decay into a ^-BPS state. 
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Now we use the BPS mass formula Eq. fl2.10|) , the formula for the decay process 
Eq. (p.2|), and and the definitions of the angles in Eq. (p^), to find after a tedious 
calculation that: 



M 4 + M 4 + M 2 4 - 2(M 2 M 1 2 + M 2 M| + M X 2 M 2 2 ) = -4r| 



sin(fl + g QP ) 

<5/?^R r- 5 m An 

sin # 



I n P„ ( mQ R f^- + nP fl ) - Sl Q fi f nF, + mQ R 

r sm v J \ sin 



2 

(3.7) 



Vanishing of the RHS is a necessary condition for marginal stability. 

This condition can be usefully rewritten by eliminating the angles 0,6qp and re- 
verting to n , r 2 coordinates for the modular parameter of the torus. It is convenient to 
introduce the following quantity depending on charges of the initial and final states as 
well as the moduli: 

E = 4= {m Sl Q\ - nr x P 2 R -{mA n)Q R ■ P R ) (3.8) 
v A 



Then we find that the equation for marginal stability is: 

mAn\ 2 ( E \ 2 1 /. x9 9 \ 

This is the "master equation" governing all two-body decays of |-BPS states in this 
theory. However we will need careful analysis to see when the equation does actually 
describe such a decay and what type of decay it describes. 

Note first of all that the equation is invariant under the transformation: 

T\ — > r 2 = —7*1, s\ — > s 2 = —Si, mi — > m 2 = m — mi, ni — * n 2 = n — rii (3.10) 

under which m An and E both change sign. This corresponds to interchange of the 
two decay products. 

If the RHS of Eq. ( |3.9| ) can be shown to be positive definite, this will be a circle in 
the torus moduli space with centre at: 

^ 11 2 ^ \ 2nsi ' 2nsi) ^ ^ 

and radius 

J(mA n) 2 + AmnrxSi + E 2 (3.12) 

2nsi 
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Because there is no restriction on the signs of r, s, it may appear that the RHS of 
Eq. ( |3.9| ) is not positive definite. However, after a little computation we are able to 
rewrite it as: 



(m A n) 2 



1 / 2 

+ Amnrisi + E 2 = —y[(m A n)Q R P R - (msi Q 2 R - nr 1 F|) cos 9 QP ] 

+(msi Q 2 R + nn Pi) 2 sin 2 9 QP ) (3.13) 



which is a sum of squares. Therefore the equation does indeed describe a nontrivial 
circle in every case. 

The next step is to check whether this circle intersects the upper half-plane. There 
are two cases. If — > then the centre of the circle is in the lower half plane. The 

si r 

circle will then intersect the upper half plane only if it intersects the real axis, which 
happens if: 

(m A n) 2 + AmnriSi > (3.14) 

It is easy to see that: 

(m A n) 2 + Amnr 1 s 1 = trF 2 - 2detF (3.15) 

where 



nvn\ nr\ 
ms\ mn\ 




(3.16) 



Now if oti, «2 are the eigenvalues of F then: 

trF 2 - 2detF = {a x - a 2 ) 2 (3.17) 

This is positive if a\, a 2 are both real, and negative if they are complex conjugate pairs. 
Therefore when is positive, only decays for which the eigenvalues of F are real can 
produce genuine curves of marginal stability in the upper half plane of r-space. 

If j- < then the circle has its centre in the upper half plane, and therefore always 
has a finite region in the upper half-plane. 
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4. Analysis of the marginal stability curves: |-BPS decay prod- 
ucts 

4.1 Equations of the curves 

To analyse the equation of marginal stability we have obtained, let us first consider the 
special case when both decay products are |-BPS. This requires that the electric and 
magnetic charge vectors of the decay products be proportional. The equation for the 
charges of the decay products Eq. Q3.1J ) can now be written: 




miQ + nP\ + / m 2 Q + r 2 P . ^ _^ 

siQ + riiP / \s 2 Q + n 2 P , 

with nii, rii, r,, Si satisfying: 

rrii + m 2 = m , n i + n 2 = n , r i + r 2 = Si + s 2 = (4.2) 

and where the electric and magnetic (upper and lower) components of each charge 
vector are proportional to each other. Note that this is the equation for the full, rather 
than projected, charge vector. The absence of any term out of the plane of Q, P comes 
from the fact that if such a term were present, it would be impossible to make the 
electric and magnetic charges proportional in both decay products. Because the above 
equation is for the full charge vectors, integrality of the charge lattice requires that 
m,,rj, Si, rii are integers. In case all four integers (for each i) have a common factor 
then the decay will be into three or more final states. Since we want to focus on 
two-body decays, we should exclude such cases. 

Proportionality of electric and magnetic charges is equivalent to requiring that the 
determinant of the associated matrices vanish: 

(let I™ 1 ri I =0 (1.3) 



and 



si ni 



. m — mi 

det ( |=° ( 4 - 4 ) 

Si 
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The first of these equations is solved by the substitution: 

(4.5) 





where a, b, c, d are defined only upto an overall reversal of sign. The second equation 
then tells us that 

mn + bcm — adn = (4-6) 

Suppose now that the original dyon was primitive, namely (m,n) = (1, 1). In this 
case Eq. ([O]) becomes 

ad - be = 1 (4.7) 

and therefore the decay products are parametrised by a matrix in PSL(2, Z). In going 
to the coefficients a, b, c, d, we see that they are invariant under the scaling a, b,c,d —>■ 
Xa, X~ 1 b, Ac, \~ l d as well as an exchange a, b,c,d —+ —b, a, —d, c. These transformations, 
along with Eq. ( |4.7| ) can be used to show that a, b, c, d are unique integers [|]. 

Making the substitutions (m, n) = (1, 1) as well as Eq. ( |4.5| ) in the curve of marginal 
stability Eq. ( |3.9| ), and using the PSL(2,Z) property, we find that the curve reduces 
to: 



ad + bc\ 2 ( E \ 2 



^ + w =7-^( 1 + ^ 2 ) ( 4 - 



led ) V 2cd J 4c 2 d 2 



where 

E = -^=(cdQ 2 + abP 2 -{ad + bc)Q-P) (4.9) 



This is the equation found by Sen in Ref, _ 

2cd au uo 2cd 



These curves are circles with centre at a ^ + , fec and radius v/ i + r f 2 ■ They intersect the 



real axis in the pair of points 

b a , . 

2-c < 410 > 

Sen showed that, for primitive dyons, two different curves never intersect in the upper 
half plane but can touch on the real axis in r-space. This implies that a given primitive 
-j-BPS dyon can at most be marginally unstable to decay into a single definite pair of 
|-BPS dyons at a given point in moduli space. 

While the fractions |, - need not in general be positive or lie between and 1, 
they can be brought into the form of positive fractions between and 1 by a modular 
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transformation. Suppose for example that | does not lie between and 1. Then for 
some suitable integer N, we define b' = b — dN such that < V < d. For the same 
N) we can show that a' = a — cN satisfies < a' < c. As a result, < — , K < 1. 
The formula for E above is unchanged under this transformation if we simultaneously 
re-define Q — > Q — NP, and the curve of marginal stability is invariant if we also send 
n -> n + AT. 

To complete the discussion of decays into |-BPS final states, we need to consider 
the case of dyons that are non-primitive, i.e. (m, n) ^ (1, 1) . In this case we can obtain 
the curve of marginal stability by starting from Eq. (|3.9|) and making the appropriate 
substitutions from Eq. ( f4.5|) and Eq. ( |4.6|) . The coefficients ad, ab, cd, be are still integers 
but they no longer describe a matrix in PSL(2, Z). Instead they satisfy the condition: 

adn — bem = ran (4-11) 

Moreover, one can check that a, b, c, d are not unique in this case. However only the 
combinations ad, ab, cd, be actually appear in the curve of marginal stability, so this 
curve is unique and can be written: 

^ _ j 2 + ( T2 + _|_y = _L_ (mV + El (4 . 12) 



E = -j= (mod Q 2 + nab P 2 - (nad + mbc)Q ■ P) (4.13) 



where 



This is the most general curve of marginal stability for decay into |-BPS dyons. 

Examining the curve we find that it intersects the real axis at the points - and 
Even though m, n are co-prime, we cannot be sure that mb, nd are co-prime, so the 
latter fraction is not necessarily reduced to lowest terms. We will discuss the geometry 
of these curves in a later subsection. 

4.2 Farey fractions and Ford circles 

In this subsection we briefly review some mathematical constructions that will facilitate 
the analysis of the |-BPS curves of marginal stability. In the mathematical literature 



one encounters the notion of a Farey sequence F n (see for example Ref. [23]). This is 



the set of all fractions (reduced to lowest terms) with denominators < n and taking 
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values in the interval between and 1, arranged in order of increasing magnitude. As 
an example we have: 

1 1 1 2 1 3 2 3 4 1, 

(4.14) 



1'5'4'3'5'2'5'3'4'5' 1 



Relevant properties of Farey sequences, for us, are the following (more details can be 
found in Ref. ||23||). For any pair of fractions 4 and - that appear consecutively in any 
Farey sequence, we have ad — be = ±1. We can always order them so that the sign is 
positive, therefore ad — bc= 1. Given any such pair, a new fraction called the mediant 
is given by: 

mediant f — , — ) = — — - (4.15) 
\d cj c + d K ' 

The mediant lies between the two members of the original pair and will occur between 

them in subsequent Farey sequences. Moreover, if we define 

h a + b 



k c + d 



(4.16) 



then hd — kb = 1 = ok — ch. Thus the fraction \ will occur after \ as well as before - 

k a c 

in some Farey sequence. 

The above construction, which is seen to be related to the structure of the discrete 
group PSL(2, Z), can be geometrically visualised in terms of circles called Ford circles. 
These will turn out to be helpful in understanding the properties of the Sen circles of 
Eq. ( |4.8|) . For a pair of co-prime integers a, c such that < - < 1, the associated 
Ford circle p3| C (|) is a circle centred at (|, ^) with radius It is tangent to the 
horizontal axis at -, and can be thought of as "sitting above" this fraction. The size of 
a Ford circle is inversely proportional to the square of the denominator of the fraction. 
Accordingly the largest possible Ford circles, above the points y and j, have radius |. 

The key property of Ford circles is that (i) two Ford circles never intersect, (ii) two 
Ford circles associated to the fractions jj and - (without loss of generality we assume the 
second fraction to be the larger one) are tangent to each other if and only if ad — bc= 1. 
In terms of Farey sequences, if two fractions are consecutive in any Farey sequence then 
they are associated to a pair of touching Ford circles. Conversely if two Ford circles 
touch then their corresponding fractions are consecutive in some Farey sequence. 

Finally we describe a construction that will be closely related to marginal decays 
of dyons. For any pair of touching Ford circles associated to fractions ^ and ^ with 
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0/1 1/4 1/3 1/2 2/3 3/4 1/1 



Figure 1: The Ford circles associated to F± 



ad — be = 1, there is another circle that (for lack of a better name) we will refer to as 
the "dual Ford circle" C(k, -) that is centred on the real axis and passes through the 
points 2 an d - on the real axis. This circle has the property that it also passes through 
the point at which the two Ford circles touch [p^j . 

4.3 Analysis of the decays: Sen circles and Ford circles 

Now let us return to the decay of a primitive |-BPS dyon into two |-BPS dyons. As 
we have seen in the previous subsection, the decay products are defined in terms of a 
matrix in PSL(2, Z). This matrix defines a pair of fractions 4 and - with ad — bc = 1. 



By the shift T\ — > T\ + N, as in the discussion below Eq. ( |4.10| ), we can make both the 
fractions lie between and 1. Now the Ford circles associated to these two fractions 
are tangent to each other. The dual Ford circle C (4, -) has its origin on the real axis 
at the midpoint of these two fractions, at ad ^ c ■ Its radius is given by half the distance 
between these two fractions, namely J^. Thus the equation of this dual Ford circle is: 



ad + be 
led 



Ac 2 d? 



(4.17) 



Comparing with Eq. 



, we see that the dual Ford circle is the limit of the Sen 



circle for marginal decays of a primitive t-BPS dyon into two ^-BPS dyons, as E — > 0. 



(Recall that E was defined in Eq. Q4.9p ). Conversely, the Sen circle can be thought of as 
a deformation of the dual Ford circle with deformation parameter E. For given integers 
a, b, c, d, both circles are centred at the same value of T\ but have their centres vertically 
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displaced from each other. The radius of the Sen circle is such that it intersects the 
real axis in the same pair of points as the dual Ford circle. Note that can be positive 
or negative, so the Sen circle can be displaced either downwards or upwards relative to 
the dual Ford circle. 

This similarity is intriguing and may point to a more profound relation between 
Sen circles and Ford circles that we have not yet uncovered (in particular, it seems 
plausible that by deforming the K3 moduli one can set E — > 0, which would make 
the two circles actually coincide). However, already the relation we have exhibited is 
sufficient to understand a key property of Sen circles, which is that they do not intersect 
in the upper half plane, but only on the real axis|§. 

This can can be seen as follows. Every Sen circle is associated to a dual Ford circle 
and thereby to a pair of Ford circles. Consider the two Sen circles associated to a, b, c, d 
and h,p, k, q with ad — be = pk — qh = 1. Clearly we have | < ^ as well as | < |. The 
two possible orderings of the fractions are |, |, ^ , | and |, f , f , |- The first ordering is 
ruled out by the Ford circle construction, since it implies that the Ford circle of the 
first fraction touches that of the third one, while the Ford circle of the second fraction 
touches that of the fourth one. This contradicts the fact that all the Ford circles are 
non-overlapping. Thus only the second ordering is possible, where we have the fractions 
i 2 , t,- in increasing order. Let us consider the case where - = 7, so that the Sen 

a ' c ' k ' q c k ' 

circles touch on the real axis. Clearly the dual Ford circles also touch on the real axis, 
which means the three fractions |, f , | are consecutive terms in a Farey sequence. 

We want to show that the Sen circles in this case do not intersect in the upper half 
plane. This imposes a condition on the slopes of the Sen circles at the real axis. From 
Eq. ( |4.8| ) we find that the slope at the real axis is given by: 

tan</> = ±-^ (4.18) 

where the two signs hold for the two intersection points. Now it is easy to check that 
the condition we are seeking is: 

E{a, b, c, d) + E(a, p, c, q) > (4.19) 

This is of course satisfied if both E's are positive, though that is not in general the case. 
But even in the general case the condition above does hold, as we now demonstrate. 
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From the definition of E one finds that: 



E(a,b,c, d)+E(a,p,c,q) = -$=(c(q + d)Q 2 + a(p + b)P 2 - (a(q + d) + c(p + b)) P ■ Q 

(4.20) 

Now we use the fact, explained in the discussion around Eq. ( |4.16| ), that if three fractions 



are consecutive in a Farey sequence then the middle one is the mediant of the other 
two. Hence we have: 

" = (4-21) 

c q+d v ' 

from which we get: 

Na = (p + b), Nc=(q + d) (4.22) 
for some integer N > 1. It follows that: 

E(a, b, c, d) + E(a, p, c, q) = —j={cQ - aPf > (4.23) 



as desired. By similar methods the non-intersecting property of Sen circles can be 
proved for the case where |, ^ , |, - are all distinct fractions. 

4.4 Analysis of the decays: non-primitive case 

The above discussion was for the case of a primitive dyon as the initial state. Now 
let us look at the case where the initial state is a non-primitive dyon. In this case the 
Sen circle is replaced by Eq. ( |4.12 ), which intersects the real axis at the points - and 
Let us now analyse the condition Eq. ( [4.1 1|) in some detail. Because m and n are 
co-prime, writing this condition as adn = m(bc + n) tells us that m divides ad and also 
that n divides be. Therefore we can rewrite Eq. (|4.11|) as: 

^ - * = 1 (4.24) 
m n 

where each of the terms on the LHS is an integer. This can only be realised if, for some 
(not necessarily prime or unique) factorisation of m and n; 

m = pq, n = kl (4-25) 

we have that: 

/ a ,/ b . c „ d * \ 

</ = -, tf= p e= v d>=- (4.26) 
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are all integers. Evidently they satisfy a'd! — b'c' = 1. Substituting in the curve of 
marginal stability for this case, Eq. ( |4.12| ), we find: 

/ p a'd' + b'c' V ( p E> V 

where 



VA \k q J 

This curve intersects the real axis at the points: 



p V pa' 
7<T ld' J 



(4.29) 



For a fixed value of ?, the set of intersection points is in one to one correspondence with 
those for the primitive case, where using Ford circles (or the methods of Ref.||) we saw 
that curves of marginal stability do not intersect. However the value of ? is not fixed. 
For given m, n specifying a non-primitive dyon, Eq. ( |4.25| ) permits several solutions for 



p and / in general. For each of them we obtain a construction in 1-1 correspondence 
with the set of curves of marginal stability for the primitive case, and it appears quite 
likely that curves from one of these sets can intersect with curves from another set. 
This would result in curves of marginal stability that intersect each other in the upper 
half plane. 

To generate examples, it is convenient to revert to the notation in which the charges 

(777*1 ^*1 \ 
J satisfying Eqs. fl4.3Q 
si «i/ 

and (|4.4|) . From these two equations we find that: 

777177, + 771777 = 77777 (4.30) 

from which we see that mi is a multiple of m. We write: 

mi = mai (4-31) 

where a± is another integer. The equations now yield the following general form for 
the matrix: 



/ 777i 


n | 


-[ 


V si 


nil 





mn qi (1— qi ) s 

mo i 



s 




(4.32) 
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(4.33) 



The strategy is now to choose a value for ai\ and then look for the set of s\ that divide 
mna,i{\ — ai). Finally to ensure that we are dealing with a two-body decay we must 
check that there is no overall common factor in either of the matrices 

m — mi —Ti 
-s\ n — rii 

In this way we can generate a large number of examples of curves of marginal stability 
for non-primitive dyons decaying into a pair of |-BPS dyons. 

To check the possible intersections of such curves, we recall that they intersect the 
real axis in the points mi ~ m and ^ . If two such intervals intersect then the curves 
will necessarily intersect in the upper half-plane. Let us consider a definite example. 
Suppose (m, n) = (2, 3). Then choosing ai = 1, we see that si can be arbitrary. On the 
other hand choosing ct\ = 2 we find that the allowed values of s\ are 1,2,3, 4, 6, 12. It 
is easy to check that for the very simplest choices the curves do not intersect. However, 
picking oli — 1, s\ — 7 and a.\ — 2, s\ — 12 we find that all the conditions are satisfied 
and the decay products are given by the matrices: 

/2 N 

ai = 1, Si 

on = 2, Sl = 12 : I * I , | ~ ' I (4.34) 

\12 -3J \-12 6J 

In terms of the integers a, b, c, d the two decay processes are parametrised by the ma- 
trices: 






(*) . = L J ( 4 - 35 ) 



Each matrix satisfies 3ad — 2bc = 6. 

Now the curves of marginal stability for the two decay modes intersect the real axis 
at the following values: 

2 

(0 n = 0, - 

(*) n = (4.36) 
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These two intervals are overlapping, hence the associated curves must intersect in the 
interior of the upper half plane. We conclude that curves of marginal stability for the 
decay of non-primitive dyons can in general intersect in the upper half plane, unlike 
what happens for primitive dyons. It would be important to understand the physical 
and mathematical reasons why the curves intersect, as well as the consequences of this 
fact. 

5. Analysis of the marginal stability curves: |-BPS decay prod- 
ucts 

5.1 Decays into a |-BPS and a ^-BPS dyon 

We now consider decays of a |-BPS dyon into one |-BPS and one ^-BPS dyon. This 
is parametrised as in Eq. ( p.2|) . If the first decay product is taken to be |-BPS then 
we must impose the condition Eq. (|4.3|) which is solved by Eq. (|4.5|) . However, the 
coefficients m;, r,, Sj, rii are no longer required to be integers and therefore nor are 
a, b, c, d. Moreover we want the second state to be |-BPS and therefore adn—bcm ^ mn. 
Finally, as indicated earlier, we must check that the curve we obtain from Eq. ( p.9| ) 
actually describes the forward and not the reverse decay process. 

Consider the case where the initial state is a primitive dyon with (m,n) = (1, 1). 
For this case we find the curve of marginal stability to be: 

( T > ~ + { T > + £)' = M (mi ~ ni? + 4riSl + £2 ) (5 ' 1) 

where 

E=^={ Sl Ql- n P 2 R - {m x - n x )Q R ■ P R ) (5.2) 
v A 

On replacing mi, r\, si, ni by their expressions in terms of a, b, c, d we can also bring it 
to the form: 

with 

E = —= (cd Q 2 R + ab P 2 R - (ad + bc)Q R ■ P R ) (5.4) 
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The equation is very similar to the Sen circle for decays of a primitive dyon into |-BPS 
decay products. However, the constraints on a, b, c, d are quite different. Instead of 
analysing this case further, we will return to it as a special case of the more general 
decay into two t-BPS states. 

5.2 Decays into two |-BPS dyons 

We now address the case in which the initial ^-BPS dyon decays into a pair of |-BPS 
dyons. Again we start with the primitive case, (m,n) = (1, 1). The relevant curve of 
marginal stability is the same as in the previous subsection, Eq. (|5.1|), except that the 



rrii Ti \ 

determinants of I ] are both nonzero. (Later we will also be able to specialise 

\ Si rii J 

to the case where one of them is zero.) 

Let us now address the constraints on the final state parameters that are required 
to ensure that the decay process corresponds to the correct branch of Eq. ( |3.6| ). First 
of all, the quantity A that appears in the BPS mass formula Eq. (|2.10| ) involves a 
square root, and we have taken all square roots to be positive. This has the following 
consequence. Observe that: 

A(m;Q + r,P, Sl Q + mP) = det 1 ) A(Q, P) (5.5) 

Positivity of A on both sides of the equation imposes the condition: 

det I 1 1 1 > 0, i = 1,2 (5.6) 



Si Tl l 



Since 



we find that: 





(5.7) 



m,\n\ — r\S\ > max (mi + n\ — 1, 0) 
For what follows, it will be convenient to introduce the eigenvalues Pi , 71 of 




1712 T2 

and /?2,72 of ( ) . Because the two matrices commute (they are of the form F 

S2 n 2l 
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and 1 — F) they can be simultaneously diagonalised, from which we see that: 



We will now examine the quantities -tj on the curve Eq. ( |5.1|) . For convenience, 
we would like to choose a particular point on the curve and evaluate these quantities 
there. The possible results are as follows. If we find ^- > 1 at a point, then the 
marginal stability curve cannot correspond to M — Mi + M 2 . It may correspond to 
either Mi = M + M 2 or M 2 = M + M\. Which of the two cases it corresponds to is 
then not very important, but can be distinguished by looking at If on the other 
hand we find -rr < 1 then we have the possibilities of being on the correct branch 
M = Mi + M 2 or on the wrong branch M 2 = M + M%. This time it is essential to 
distinguish the two, which can again be done by evaluating Being on the correct 
branch requires ^ < 1 for both i — 1 and 2. 

In any of these cases, having determined the relevant branch of Eq. ( |3.6|) at one 
point on the curve, we can be sure that we will not cross over to another branch 
elsewhere on the same curve, since crossing from one branch to another requires passing 
through a point where one of the masses vanishes. But the BPS mass formula does not 
vanish for any value of the moduli, so this is not possible (unless the charges of that 
state vanish identically). 



the real axis. This will happen if the eigenvalues fli,ji are both real (without loss of 
generality we take 71 > (3i). Then, a convenient point at which to evaluate the mass 
ratios is one of the intersection points of the curve with the real axis. Setting r 2 = in 
Eq. ( |3.9| ), we get the following equation for t\\ 



Pi + fa = 1 = 7i + 72 



(5.9) 



From the determinant conditions above, we have: 



#l7i>0, (1-A)(1-7i) >0 



(5.10) 




n 

ni = rrii — T iSi 

Tl 



(5.11) 



M 2 

Now let us consider the expression -j-k in the limit r 2 — > 0. We have: 



M{ 
M2 



[{mxQ R + t x Pr) - n(siQ R + uxPr)} 



2 



r 2 ^0 [Qr - T\P R \ 2 

n \ r. 
R 



[K-nsOgfl-riK-a^ia 



[<3* - t-iPr] 

Using Eq. ( 5 . 11|) we now get: 



(5.12) 



Mi 
~M 



I mi — TiSi| (5.13) 



On the real axis, Eq. (|5.1|) gives: 

n = ^- ( ± (71 - A) I + (mi - m)) (5.14) 

Inserting this into Eq. (|5.13|) we find: 

mi — T\Si = j3\ or 71 (5.15) 

Let us first consider the case mini — ri«i > 1. We will show that in this region the 
decay is not the desired one, but corresponds instead to a branch of Eq. (|3.6| ) describing 
a reverse decay. With this condition on the determinant, one of the eigenvalues (say 
71) must be > 1. Positivity of the second determinant, which equals (1 — #l)(1 — 71), 
tells us that if 71 > 1 then also (3\ > 1. Thus we have that both eigenvalues are > 1. 
It follows that ^ > 1 and we are, as promised, on the wrong branch. 

Next suppose — riSi = 1. The above considerations then show that f3\ = 

71 = 1. Then we ^ = 1. This means M 2 = and therefore the charges associated to 

(mi r x \ [0 o\ 

the second state are identically zero. In other words, = • This is a 

\ai ni) \0 Oj 

trivial case where the first decay product is the original state itself. 

Let us note at this point that if mi,ri,Si,ni had been taken to be integers, and 
the corresponding state was restricted not to be |-BPS, we would necessarily have 
mini —T\S\ > 1. We have shown that all such cases do not correspond to a valid decay 
of M into Mi and M2, therefore no such decays exist for integer coefficients. This is 



one of the key results of Ref.[14 . 
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That leaves the case 



< m\ri\ — T\S\ < 1, < m 2 n 2 — r 2 s 2 < 1 (5.16) 

which can only be satisfied for fractional coefficients. 

Requiring j3\^\ < 1 and also P2I2 = (1 — — 7i) < 1 we see that < A, 71 < 1 
and < (3 2 , 72 < 1- From this and Eq. ( |5.15| ) we find ^ < 1, ^ < 1 and this indeed 
corresponds to the decay process that we were looking for. Thus Eq. (|5.16|) provides a 
necessary condition for the coefficients m-i, n%, r\, s\ in Eq. ( |3.2[) in order to have a decay 
of the original dyon into two |-BPS dyons. Under this condition, our curve Eq. ( |3.9|) 
describes the marginal stability locus in the Ti, r 2 plane. However this is a locus of co- 
dimension 2 in the full moduli space, for the following reason. Fractional mi,ri, S\,rii 
means that the decay process in terms of the original integral charge vectors was into 
states living outside the Q, P plane. This is precisely the case, referred to earlier, 
where the moduli in M need to be adjusted to make the final state charges (after R 
projection) lie in the same plane as the initial state charges It remains to find a 
sufficient condition on the values of mi, ri, s\, n\ as well as to understand more precisely 
the condition on the moduli matrix M which put the projected charge vectors in the 
plane of the decaying dyon. 



6. Discussion 

We have found a general equation for marginal stability of |-BPS dyons to decay into 
two final state particles, Eq. ( |3.9| ). Analysis of the equation reveals many distinct cases 
with different properties. We believe this analysis can be easily extended to multi- 
particle final states. The construction of Ford circles and especially their dual circles 
proved useful in this analysis and we suspect that there may be a deeper mathematical 
relationship to the Sen circles of marginal instability for primitive dyons decaying into 
|-BPS final states. 

To complete the discussion of the previous section of decays into j-BPS states we 

(mi r\ \ 
J are complex, in which 
Sl Tli J 

case the curve of marginal stability does not intersect the real axis. Also we need to 
generalise to the case of non-primitive initial states. 
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We see that decays into |-BPS final states are labelled by integers a, b, c, d that 
are in PSL{2, Z) for primitive initial states and obey the more complicated relation 
Eq. ( |4.6| ) when the initial state is a non-primitive dyon. However, such integers do not 
occur when the final state consists of |-BPS dyons. It would be nice to understand the 
physical origin of this PSL(2, Z) and its generalisations in the cases where they occur. 



We expect this can be done through the string network|20], [21], ^2| representation of 
the dyons. 

The role of non-primitive dyons and their decays has been somewhat mysterious 3 
since the observation of Ref . that the quantum states of such dyons are not counted 



by the famous genus-2 modular form of Ref. ||T3| but appear to be connected to a higher- 
genus Riemann surface. We have exhibited the |-BPS decays of such dyons and noted 
that their curves of marginal stability intersect in the interior of moduli space. This 
may be helpful in resolving the puzzle of their role in the counting problem. 
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